Open Access Journal of Sports Medicine 



Q Open 



Access Full Text Article 



Dovepress 

open access to scientific and medical research 

ORIGINAL RESEARCH 



Mathematical simulation of energy expenditure 
and recovery during sprint cross-country skiing 



This article was published in the following Dove Press journal: 
Open Access Journal of Sports Medicine 
14 May 2014 

Number of times this article has been viewed 



John F Moxnes 1 
Eldbjorg Dirdal Moxnes 2 

'Protection and Societal Security 
Division, Norwegian Defence 
Research Establishment, Kjeller, 
Norway; 2 Department of Geosciences, 
University of Oslo, Oslo, Norway 



Correspondence: Eldbjorg Dirdal Moxnes 
Department of Geosciences, 
University of Oslo, PO Box 1 072 
Blindern, Oslo, Norway 
Email eldbjorg_moxnes@hotmail.com 

John F Moxnes 

Protection and Societal Security 
Division, Norwegian Defence Research 
Establishment, PO Box 25, 
2007 Kjeller, Norway 
Email john-f.moxnes@ffi.no 



Purpose: A cross-country sprint competition relies on maximal effort durations of 3-4 minutes. 
Significant anaerobic energy contribution is expected. Anaerobic energy contribution has been 
estimated in different sports to date from the accumulated 0 2 deficit. However, the 0 2 -deficit 
model can be questioned. We investigate anaerobic energy contribution by applying other 
methods than the 0 2 deficit. 
Methods: Theoretical model development. 

Results: For sprint cross-country competitions, the anaerobic energy contribution was 20%-25% 
independent of the employed mathematical model. Recovery times of a minimum 20 minutes 
were found to be required after sprint races to be sure that the performance in subsequent heats 
was not influenced. 

Conclusion: The 0 2 -deficit model gave anaerobic energy results in agreement with other models 
from the literature. Recovery times of a minimum 20 minutes were found to be required after 
sprint races to be sure that the performance in subsequent heats was not influenced. 
Keywords: aerobic, anaerobic, models, endurance sport, recovery 

Introduction 

In endurance sports, energy expenditure and the ability to utilize metabolic energy 
to produce external work are the two main performance-determining factors. Energy 
expenditure is dependent on the body's ability to synthesize (produce) and con- 
sume adenosine triphosphate (ATP). Maximal use of ATP in mammalian muscles is 
around 1.7xl0~ 5 mol ATP/g tissue/second. ATP stores in skeletal muscle tissue are 
around 5x1 0~ 6 mol/g tissue. Therefore, maximal use can only be applied for around 
0.5 seconds without production of ATP. 1 For a muscle mass of, say, 10 kg, this gives 
170 mmol/second, which is much more than around 1 mmol/second during rest. 2 
Despite the more than 100-fold increase in the rate of ATP consumption from rest to 
maximal-intensity exercise, the energetic demands of muscles in endurance sports are 
usually satisfied without depleting intracellular ATP. 3 Three sources are available for 
ATP production: 1 ) ATP can be produced aerobically in the mitochondria by oxidative 
phosphorylation (aerobic energy), 2) ATP can be produced by anaerobic synthesis due 
to glycolysis/glycogenolysis (lactic energy), and 3) ATP can be produced by phos- 
phocreatine (PCr) breakdown to creatine (Cr) (adenosine diphosphate + PCr gives 
ATP + Cr in the creatine kinase [CK] reaction) (alactic energy). The maximal amount 
of anaerobic energy that can be utilized is proportional to the sum of the maximum 
amount of Cr and lactate that can be accumulated in the body 1 The equilibrium con- 
stant of the CK reaction is around 20. 3 Therefore, the slightest drop in ATP allows the 



submit your manuscript | www.dovepi ess.com 
Dovepress 

http://dx.doi.org/ 1 0.2 1 47/OAJSM.S62020 



Open Access Journal of Sports Medicine 20 1 4:5 115-121 



115 



©2014 Moxnes and Moxnes This work is published by Dove Medical Press Limited, and licensed under Creative Commons Attribution Non Commercial (unported. v3 0) 
| License. Ihe full terms of the License are available at httpy/creativecommons.org/licenses/by-nc/i.O/ Non-commercial uses of the work are permitted without any further 
permission from Dove Medical Press Limited, provided the work is properly attributed. Permissions beyond the scope of the License are administered by Dove Medical Press Limited. Information on 
how to request permission may be found at http //www dovepress com/permissions php 



Moxnes and Moxnes 



Dovepress 



reaction to proceed to ATP, and the ATP concentration stays 
nearly constant until almost all the PCr is exhausted. 3 

The aerobic energy release is the main source in endur- 
ance sports with durations above 5 minutes. Anaerobic 
energy release is the main source with competition times 
below 1 minute. 1 However, the energy release in so-called 
middle-distance sports (defined as racing times between 
1 and 5 minutes) is determined by a more complex interaction 
between the various energy sources. When such competition 
times are performed multiple times within the same day, as 
in sprint cross-country skiing or some of the speed-skating 
disciplines, a rapid recovery of the energy systems is required. 4 
Although several studies have examined many of these aspects 
in an isolated fashion through various experiments, mathemati- 
cal modeling may provide a more integrated understanding. 

A crucial metabolic pathway contributing to muscle- ATP 
regulation in middle-distance sports is glycolysis/glycogenol- 
ysis. The end point of glycolysis is pyruvate, which represents 
a metabolite that can be reduced to form lactate or oxidized to 
C0 2 or H,0. The blood lactate concentration is the result of 
the appearance and disappearance of lactate. 5 - 6 The working 
muscles and various tissues produce lactate and release it 
into the plasma. Both resting and submaximal working of the 
skeletal muscles, as well as the brain, heart, liver, and kidney, 
remove lactate from circulation by the use of mitochondria. 7 " 10 
It has been suggested that lactate acts as an intermediary 
for the transport of carbohydrate from cells and tissue with 
relatively low oxidative capacity, through the blood lactate 
pool to cells and tissues with high oxidative capacity. 7-10 
Oxidation in the exercising skeletal muscles accounts for 
around 70%-80%, and the liver (gluconeogenesis in the 
liver by the Cori cycle) accounts for around 20%-30% of 
the lactate disappearance during submaximal exercise. 11 
Around 20% of the carbohydrate oxidized during exercise 
passes through the blood lactate pool before being oxidized to 
carbon dioxide. 11 During mild or moderate exercise (typically 
below peak fat consumption around 50%-60% of maximal 
rate of oxygen consumption [V0 2ma J), the rate of glycolysis 
increases several-fold, but only leads to a moderate lactate 
accumulation in the body. 1 During exercise, >60% of V0 2max 
(around peak fat consumption) 2 up to the lactate threshold 
(LT), 1 (above the LT, no steady state of lactate is developed 
through time), the steady-state lactate concentration reaches 
significantly higher values than at rest. 1 Middle-distance 
sports are performed above the LT, and a progressive accu- 
mulation of lactate in muscles and blood develops. In such 
cases, the capillary blood pH can fall from 7.45 to around 
7. 05, 12 whereas for skeletal muscles the pH can fall from 7.1 



to 6.4. 13 It is not lactate itself that leads to muscle fatigue, 
but most likely the associated fall in pH. 14,15 

After step changes in work rates for moderate-intensity 
exercise below the LT, PCr levels, as aerobic power, develop 
exponentially and reach a steady-state level. Below the LT, 
a strong similarity has been reported for the time constant 
for 0 2 kinetics response and PCr consumption. 16 Above 
the LT the anaerobic glycolytic energy supply becomes 
significant, and the similarities between the rate constants 
have not been systematically reported. 17 PCr recovery is 
mainly due to oxidative ATP synthesis. However, even when 
oxygen is excluded, PCr stores may be rebuilt by anaerobic 
glycolysis. 18-20 In general PCr is commonly acknowledged as 
an energy buffer, which supports the transient failure of other 
metabolic pathways to support ATP. During recovery of PCr 
levels, the lactate (the pH must be returned) and adenosine 
diphosphate must be returned to a normal state. 

In this paper, the different anaerobic energy contributions 
by applying different models are studied. It is found how 
the anaerobic portion of total energy depends on time, and 
a sprint cross-country skiing scenario is used to study how 
racing time and recovery depends on the concentration of 
lactate in the body. Cross-country skiing was chosen because 
it is a whole-body exercise where both anaerobic and aerobic 
energy expenditure is important. 

Sprint cross-country skiing is a relatively new racing form 
that demands high aerobic and anaerobic power along with 
highly developed technical and tactical skills. 21 The heats 
in sprint skiing (eg, two different semifinals) are performed 
in sequences after each other, which may lead to individual 
differences in the duration of breaks before the final. If the 
duration of breaks between heats is too short, this rationally 
has an impact on the final sprint performance. This topic was 
discussed after the Vancouver Olympics in 2010, where the 
duration of breaks from the two semifinals to the final was 
around 10 and 15 minutes, respectively. 

The main hypotheses of the present study were: 

1 . The estimated anaerobic energy contribution during sprint 
racing is dependent on the mathematical models in the 
literature. 

2. Recovery times of a minimum 1 5-20 minutes are required 
after cross-country sprint races to be sure that perfor- 
mance in subsequent heats is not influenced. 

Methods 

Theoretical model development 

E is the aerobic energy used per unit of mass (that means 
where the used ATP is synthesized aerobically), E is the 
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anaerobic energy per unit of mass due to glycogenolysis or 
glycolysis, and E CK is the anaerobic energy per unit of mass 
due to the CK reaction. Therefore: 



E. + 



(1) 



total 
energy 



aerobic 
energy 



anaerobic energy, 
wbere ATP synthesized 
from 



glycogenolysis/glycolysis CK 



anaerobic 
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where mod means a testable model. During the race we assume 
that the two anaerobic capacities are fully exploited, since the 
racing time is above 1 minute. 22 Equation 1 can be written as: 



^ A 

E 



1- 



(2) 



Therefore, to calculate the portion of aerobic energy, 
knowledge of the lactic energy (E G ) and alactic energy (E CK ) 
is needed. 

Lactic energy 

The relationship between oxygen consumption and blood 
lactate levels (ie, the lactic energy) can be calculated accord- 
ing to di Prampero and Ferretti 23 and used to estimate the 
anaerobic energy from glycolysis only. The mL/kg 0 2 that 
corresponds to 1 mmol/L varies from 2.7 to 3.3, due to varia- 
tion in the relative blood volume versus active muscle mass 
for an athlete. 23 If this value is set to an average of 3.0 mL/kg 
0 2 , a change in the blood lactate concentration from 0.5 to 
15 mmol/L gives: 



E = 14.5x3.0 = 43.5 mL/kg 



(3) 



Alactic energy 

Alactic energy is more complicated and different models are 
in use in the literature. One model simply says that alactic 
power - Q CK (t) - is proportional to the rate of change of 
aerobic power - Q A (t) 17 > 24 - 25 - to read: 



mod 

Q CK (t) = 0Q A (t ) 



(4) 



where 9 is a constant of the proportionality parameter that 
follows. By time integration, Equation 4 gives: 



Model l:E a =e(Q A (i)-QJt^ 



(5) 



Here, the power due to the rate of 0 2 consumption (V0 2 ) is 
subtracted at the end and start of the exercise, and the number 
is multiplied by 0. It can be argued that: 25 



where rj A rj CK is the efficiency when producing ATP 
aerobically and from the CK reaction, respectively, and r A 
is a time parameter quantifying the time before the aerobic 
power reaches a steady state during a steady-state work 
rate. Typical values for r A are 30-36 seconds for moderate 
intensity exercise. 6,22 - 23,26 ~ 30 For example, di Prampero 22 sug- 
gested that r A equals 10-24 seconds. Cerretelli et al 29 found 
that % A increases linearly with concentration of lactate up to 
36 seconds, and Binzoni et al 30 found 23 seconds for all work 
rates. We set r A to 30 seconds as a compromise. 

Another method to calculate the energy due to the CK 
reaction is a variant of the so called 0 2 -deficit assumption, 25 
here called Model 2. It is assumed that for a given type of 
exercise, the rate of ATP consumption is the same for the 
same work rate - P. 

Model 2: tjtj a E a + Tffl CK E CK + tjtj 0 E g = Tjij^ 
( \ 



. 0 2 deficit j 



(7) 



Here, r\ is the efficiency during muscle contraction, 
while r\ G is the efficiency when producing ATP by anaerobic 
glycolysis/glycogenolysis. For exercise intensities exceeding 
the maximal rate of 0 2 consumption, the virtual steady state 
is the steady state that would be attained if it was possible 
to carry out the exercise under purely aerobic conditions. 23 
Obviously, this tentative (or virtual) steady state rate Evir IS 
never reached as the increase in 0 2 uptake approaches the 
maximal rate of 0 2 consumption. Here, E vir (t ) = J E vir (u)du . 
It is assumed that for a given type of exercise, J? =a + bP, 
where P is the work rate. The parameters a and b are deter- 
mined by fitting the steady-state rate of 0 2 consumption to 
work rate P. The 0 2 -deficit model applies tj a = tj ck = r/ G . 
Therefore, E CK = (E vw = E A )-E a . However, ?j A = 7j CK = t] G 
seems to have no justification in the literature. Now, we set 
rj A = fJ G = 0.6,// CK = 0.95 as the baseline. 25 

If E vtr - E A = 76 mL/kg, this gives: 31 

0,-deficitmethod:iw = — x 76 - — x 43.5 
2 CK 0.6 0.6 

= 32.5 mL/kg 
Model 1: E CK = 6>{Q A (t)-Q A (f 0 )) = 24.9 mL/kg 

Model 2:^ =^x 76-^x43. 5 = 20.5 mL/kg (8) 
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(6) 



Obviously, the results will be dependent on the 
efficiencies. 
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For comparison, E CK numbers from the 0,-deficit method 
and the two models can be compared to the number calculated 
when using the estimate by di Prampero. 22 For exercise lead- 
ing to exhaustion from 50 seconds to 10 minutes, the sum of 
the two anaerobic energy sources is essentially constant and 
equal to 1 ,400 J/kg =70 mL 0 2 /kg if it is assumed that 1 mL 
0 2 corresponds to 20 J. Using Equation 3 and subtracting 
E a =43.5 mL/kg gives E CK =26.5 mL/kg. This number is in 
agreement with the numbers in Equation 8. 

The results can be summarized as: 



Ea _ y E G + E ck 



Model 2: E a +E CK 



0 2 -deficit method: E G + E CK = E vir - E A 
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-(e»-e a ) 



(9) 



Results 

A skier utilizes the average rate of aerobic power of Q A . 
This gives: 



Ej_ _ ^ E a + E CK 



QJ 



QJ + E G +E C 



(10) 



where t is the time of the race. Ninety percent of the peak 
rate of oxygen consumption measured in the laboratory is 



used in the calculations as average aerobic power during a 
sprint race. This assumption is based on studies showing 
peak values of around 95% of their maximal rate of oxygen 
consumption, 32 as well as mean values of around 95% of the 
peak values during a sprint race. 21 

Figure 1 shows the solution of Equation 10 when Q 4 
=0.9 x 77.8 mL/kg/minute =70 mL/kg/minute. The mean 
duration of a cross-country sprint ski heat is 3 minutes 
30 seconds, and 3 minutes 10 seconds for men and women, 
respectively. 4 

An overall mean of 3 minutes 20 seconds is used further. 
Figure 1 shows the portion of aerobic and anaerobic energy 
of total energy as a function of time, and in the specific case 
of 3-minute 20-second duration, the anaerobic contribution 
is around 20%-25%. 

PCr levels quickly recover after a race (after around 30 
seconds 3 ). However, recovery from anaerobic glycolysis 
needs longer time periods. Assuming the values of the total 
energy in a sprint heat for a world-class sprint skier after 
breaks of different duration following a sprint heat, blood 
lactate level is calculated. The mean aerobic power was 
regarded as 70 mL/kg/minute, and the mean peak blood 
lactate level 15 mmol/L. After 6-, 10-, 15-, and 20-minute 
breaks after the sprint, the skier would start with 9.5, 6.5, 
3.5, and 1.8 mmol/L lactate, respectively. 32 Equation 10 
can be reformulated to study the time for a race as a func- 
tion of the initial blood lactate concentration of the skier. 
During a race of around 3 minutes 20 seconds, the energy 
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Figure I Portion of aerobic energy to total energy as a function of time (t) in seconds (s). Q A - 70 mL/kg/minute. 
Notes: Red, E G + E CK =76 mL/kg; blue, £ G + E CK -68.4 mL/kg; green, £ G + E CK -64.0 ml_/kg. 
Abbreviations: £ G , lactic energy; £ CK , atactic energy. 
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used is called E„ where denotes sprint skiing. The energy 
is then: 



QJ 



E s Qa 1 + Eq + E ck 
E s = Q A t + E G + E CK 

= 70/60 x (3x60 + 20)+ 68.4 = 302 mL/kg (11) 

Let C denote the blood lactate concentration in millimoles 
per liter at the start of the race, and let / = t(C) be the time 
for the race. It can then be written that: 

t(C) = {E s -E a (C)-E CK )lQ A 
E G (C) = (15 mmol/L -C)xA, 

A = 3 mL/kg /(nmol/L) (12) 

This gives the time relative to the time when the initial 
concentration is 1.8 mmol/L (corresponding to 20 minutes' 
recovery) as: 

t(C) _ E s -E a (C)-E CK 



r(l.8) E s -E a {L8)-E CK 

, ,C- 1.8 mmol/L 
= l + A 



(13) 



As an example, using the approximation above and 
an initial concentration of 6.5 mmol/L (corresponding to 
10 minutes' recovery) gives f(6.5)/<(1.8) = 1 + 3(6.5 - 1.8)/ 
302 = 1.05. Therefore, the difference between 20 minutes' 



break and 10 minutes' break is around 5%. The second 
term of Equation 1 3 can be plotted to study the effect more 
carefully. Figure 2 shows the result. The exact result was 
found to be closer to 7.5%. 

The current findings reveal differences of around 6%, 3%, 
and 2% in racing time between breaks of duration 10 versus 
20 minutes, 10 versus 15 minutes, and 15 versus 20 minutes, 
respectively. This means that the duration of breaks may 
impact performance in the subsequent heat when the breaks 
are below 20 minutes between the heats. With 20-minute 
breaks between heats, world-class skiers show blood lactate 
levels around 2 mmol/L, which is closer to the resting levels 
of blood lactate concentration - around 0.5 mmol/L. However, 
these effects only concern low-altitude competitions. 

Conclusion 

Anaerobic energy in sprint cross-country skiing was inves- 
tigated by using two mathematical models that have been 
recently published in the literature. The models were used to 
study the anaerobic portion of energy as a function of time. 
These models have been used to study sprint-skiing racing 
scenarios in more detail. It was shown that the anaerobic 
versus aerobic portion of energy should be considered when 
choosing the duration of a sprint race; today's races contain 
between 25% and 30% anaerobic and 75%-80% aerobic 
energy. Also, the duration of breaks between different racing 
heats in sprint skiing is of importance: based on literature 
values for blood lactate during recovery, the simulations 
showed that recovery times of 20 minutes between heats of 
world-class skiers are necessary to ensure that the breaks do 
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not impact performance in subsequent heats. Only hypothesis 
2 was accepted. 

The simulations showed recovery times of 20 minutes 
between heats of world-class skiers are necessary to ensure 
that the breaks do not impact performance in subsequent 
heats. At higher altitudes, such as in the Olympic Games in 
Sochi in 2014, at 1,500 m above sea level, the resting time 
probably requires even longer durations. 

Only mean values of energy during the last 3 minutes 
20 seconds of a race were calculated, despite the possible 
variable work rates and intensities during the sprint heats. 
How the latter aspects of sprint skiing, as well as the dif- 
ferences in tactics within the heats, may impact the results 
were not considered. 

It is also assumed that aerobic power remains constant 
across differences in pH (here indicated by blood lactate 
concentration). This assumption is not obviously correct, 
because beta oxidation in the mitochondria is shown to 
be affected by pH. Earlier studies of sprint skiing showed 
no significant differences in mean or peak oxygen uptake 
across heats. 2133,34 Finally, other factors affect recovery, 
such as muscular and neuromuscular fatigue, but were not 
considered in these analyses. Overall, this paper provides 
novel knowledge about energy demands and recovery in 
cross-country sprint skiing. 
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